Importing the data-set
# National Parks in California
ca = read_csv("https://raw.githubusercontent.com/ScienceParkStudyGroup/r-lesson-based-on-ohi-data-training/gh-pages/data/ca.csv")
Rows: 789 Columns: 7
── Column specification ───────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (5): region, state, code, park_name, type
dbl (2): visitors, year
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.
head(ca)
# Acadia National Park
acadia <- read_csv("https://raw.githubusercontent.com/ScienceParkStudyGroup/r-lesson-based-on-ohi-data-training/gh-pages/data/acadia.csv")
Rows: 98 Columns: 7
── Column specification ───────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (5): region, state, code, park_name, type
dbl (2): visitors, year
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.
head(acadia)
# Southeast US National Parks
se <- read_csv("https://raw.githubusercontent.com/ScienceParkStudyGroup/r-lesson-based-on-ohi-data-training/gh-pages/data/se.csv")
Rows: 453 Columns: 7
── Column specification ───────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (5): region, state, code, park_name, type
dbl (2): visitors, year
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.
head(se)
# 2016 Visitation for all Pacific West National Parks
visit_16 <- read_csv("https://raw.githubusercontent.com/ScienceParkStudyGroup/r-lesson-based-on-ohi-data-training/gh-pages/data/visit_16.csv")
Rows: 17 Columns: 7
── Column specification ───────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (5): region, state, code, park_name, type
dbl (2): visitors, year
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.
head(visit_16)
# All Nationally designated sites in Massachusetts
mass <- read_csv("https://raw.githubusercontent.com/ScienceParkStudyGroup/r-lesson-based-on-ohi-data-training/gh-pages/data/mass.csv")
Rows: 13 Columns: 7
── Column specification ───────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (5): region, state, code, park_name, type
dbl (2): visitors, year
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.
head(mass)
Plot-2 –> ggplot
Scatter plot
scatter_plot = plot_ly(data=ca, x=~year, y=~visitors,color = ~park_name, type='scatter',mode='markers') %>%
layout(
title= list(text = "<b>Body weight vs Brain weight"),
legend = list(title = list(text ='<b>Animals')),
xaxis = list(title = list(text ='<b>Brain Weight')),
yaxis = list(title = list(text ='<b>Body Weight')))
scatter_plot
Warning in RColorBrewer::brewer.pal(N, "Set2") :
n too large, allowed maximum for palette Set2 is 8
Returning the palette you asked for with that many colors
Warning in RColorBrewer::brewer.pal(N, "Set2") :
n too large, allowed maximum for palette Set2 is 8
Returning the palette you asked for with that many colors
Warning in RColorBrewer::brewer.pal(N, "Set2") :
n too large, allowed maximum for palette Set2 is 8
Returning the palette you asked for with that many colors
Warning in RColorBrewer::brewer.pal(N, "Set2") :
n too large, allowed maximum for palette Set2 is 8
Returning the palette you asked for with that many colors
Barchart
r_group_barchart = data.table::melt(ca, id.vars='park_name') %>%
plot_ly(x = ~park_name, y = ~value, type = 'bar', name = ~variable, color = ~variable) %>%
layout(
title= list(text = "<b>Total Distribution based on Vore"),
legend = list(title = list(text= '<b>Aniamal Feature')),
xaxis = list(title = list(text ='<b>Vores')),
yaxis = list(title='Count', text='<b>Count'), barmode = 'group')
Warning in data.table::melt(ca, id.vars = "park_name") :
The melt generic in data.table has been passed a spec_tbl_df and will attempt to redirect to the relevant reshape2 method; please note that reshape2 is deprecated, and this redirection is now deprecated as well. To continue using melt methods from reshape2 while both libraries are attached, e.g. melt.list, you can prepend the namespace like reshape2::melt(ca). In the next version, this warning will become an error.
r_group_barchart
pie-chart
df_order = data.frame(table(ca$park_name))
df_order
pie_chart = plot_ly(type='pie', labels=df_order$Var1, values=df_order$Freq,
textinfo='label+percent',insidetextorientation='radial') %>%
layout(
title= list(text = "<b>Order Distributions"),
legend = list(title = list(text= '<b>Order')))
pie_chart
histogram_plot = plot_ly(data = ca, x = ~(log(visitors)), name=~code,type="histogram") %>%
layout(
title= list(text = "<b>Total sleep time of Animals based on Vore"),
legend = list(title = list(text= '<b>Vore')),
xaxis = list(title = list(text ='<b>Visitors')),
yaxis = list(title = list(text ='<b>Count')))
histogram_plot
Donut Chart / Open Pie-Chart
df_vore = data.frame(table(ca$code))
df_vore
donut_chart = plot_ly(labels=df_vore$Var1, values=df_vore$Freq,
textinfo='label+percent') %>%
add_pie(hole = 0.6) %>%
layout(
title= list(text = "<b>Order Distributions"),
legend = list(title = list(text= '<b>Order')))
donut_chart
LS0tCnRpdGxlOiAiUiBOb3RlYm9vayIKb3V0cHV0OiBodG1sX25vdGVib29rCi0tLQoKIyMgSW1wb3J0aW5nIHRoZSBsaWJyYXJpZXMKYGBge3J9CmxpYnJhcnkodGlkeXZlcnNlKQpsaWJyYXJ5KHBsb3RseSkKYGBgCgojIyBJbXBvcnRpbmcgdGhlIGRhdGEtc2V0CmBgYHtyfQojIE5hdGlvbmFsIFBhcmtzIGluIENhbGlmb3JuaWEKY2EgPSByZWFkX2NzdigiaHR0cHM6Ly9yYXcuZ2l0aHVidXNlcmNvbnRlbnQuY29tL1NjaWVuY2VQYXJrU3R1ZHlHcm91cC9yLWxlc3Nvbi1iYXNlZC1vbi1vaGktZGF0YS10cmFpbmluZy9naC1wYWdlcy9kYXRhL2NhLmNzdiIpIAoKaGVhZChjYSkKYGBgCgpgYGB7cn0KIyBBY2FkaWEgTmF0aW9uYWwgUGFyawphY2FkaWEgPC0gcmVhZF9jc3YoImh0dHBzOi8vcmF3LmdpdGh1YnVzZXJjb250ZW50LmNvbS9TY2llbmNlUGFya1N0dWR5R3JvdXAvci1sZXNzb24tYmFzZWQtb24tb2hpLWRhdGEtdHJhaW5pbmcvZ2gtcGFnZXMvZGF0YS9hY2FkaWEuY3N2IikKCmhlYWQoYWNhZGlhKQpgYGAKCmBgYHtyfQoKIyBTb3V0aGVhc3QgVVMgTmF0aW9uYWwgUGFya3MKc2UgPC0gcmVhZF9jc3YoImh0dHBzOi8vcmF3LmdpdGh1YnVzZXJjb250ZW50LmNvbS9TY2llbmNlUGFya1N0dWR5R3JvdXAvci1sZXNzb24tYmFzZWQtb24tb2hpLWRhdGEtdHJhaW5pbmcvZ2gtcGFnZXMvZGF0YS9zZS5jc3YiKQoKaGVhZChzZSkKYGBgCgpgYGB7cn0KIyAyMDE2IFZpc2l0YXRpb24gZm9yIGFsbCBQYWNpZmljIFdlc3QgTmF0aW9uYWwgUGFya3MKdmlzaXRfMTYgPC0gcmVhZF9jc3YoImh0dHBzOi8vcmF3LmdpdGh1YnVzZXJjb250ZW50LmNvbS9TY2llbmNlUGFya1N0dWR5R3JvdXAvci1sZXNzb24tYmFzZWQtb24tb2hpLWRhdGEtdHJhaW5pbmcvZ2gtcGFnZXMvZGF0YS92aXNpdF8xNi5jc3YiKQoKaGVhZCh2aXNpdF8xNikKYGBgCgpgYGB7cn0KIyBBbGwgTmF0aW9uYWxseSBkZXNpZ25hdGVkIHNpdGVzIGluIE1hc3NhY2h1c2V0dHMKbWFzcyA8LSByZWFkX2NzdigiaHR0cHM6Ly9yYXcuZ2l0aHVidXNlcmNvbnRlbnQuY29tL1NjaWVuY2VQYXJrU3R1ZHlHcm91cC9yLWxlc3Nvbi1iYXNlZC1vbi1vaGktZGF0YS10cmFpbmluZy9naC1wYWdlcy9kYXRhL21hc3MuY3N2IikKCmhlYWQobWFzcykKYGBgCgojIyBQbG90LTIgLS0+IGdncGxvdAojIyMgU2NhdHRlciBwbG90CmBgYHtyfQpzY2F0dGVyX3Bsb3QgPSBwbG90X2x5KGRhdGE9Y2EsIHg9fnllYXIsIHk9fnZpc2l0b3JzLGNvbG9yID0gfnBhcmtfbmFtZSwgdHlwZT0nc2NhdHRlcicsbW9kZT0nbWFya2VycycpICU+JQogICAgICAgICAgICAgICAgICAgICAgIGxheW91dCgKICAgICAgICAgICAgICAgICAgICAgICAgICB0aXRsZT0gbGlzdCh0ZXh0ID0gIjxiPkJvZHkgd2VpZ2h0IHZzIEJyYWluIHdlaWdodCIpLAogICAgICAgICAgICAgICAgICAgICAgICAgIGxlZ2VuZCA9IGxpc3QodGl0bGUgPSBsaXN0KHRleHQgPSc8Yj5BbmltYWxzJykpLCAKICAgICAgICAgICAgICAgICAgICAgICAgICB4YXhpcyAgPSBsaXN0KHRpdGxlID0gbGlzdCh0ZXh0ID0nPGI+QnJhaW4gV2VpZ2h0JykpLAogICAgICAgICAgICAgICAgICAgICAgICAgIHlheGlzICA9IGxpc3QodGl0bGUgPSBsaXN0KHRleHQgPSc8Yj5Cb2R5IFdlaWdodCcpKSkKc2NhdHRlcl9wbG90CmBgYAoKIyMjIEJhcmNoYXJ0CmBgYHtyfQpyX2dyb3VwX2JhcmNoYXJ0ID0gZGF0YS50YWJsZTo6bWVsdChjYSwgaWQudmFycz0ncGFya19uYW1lJykgJT4lCnBsb3RfbHkoeCA9IH5wYXJrX25hbWUsIHkgPSB+dmFsdWUsIHR5cGUgPSAnYmFyJywgbmFtZSA9IH52YXJpYWJsZSwgY29sb3IgPSB+dmFyaWFibGUpICU+JQogICAgICBsYXlvdXQoCiAgICAgICAgICB0aXRsZT0gbGlzdCh0ZXh0ID0gIjxiPlRvdGFsIERpc3RyaWJ1dGlvbiBiYXNlZCBvbiBWb3JlIiksCiAgICAgICAgICBsZWdlbmQgPSBsaXN0KHRpdGxlID0gbGlzdCh0ZXh0PSAnPGI+QW5pYW1hbCBGZWF0dXJlJykpLAogICAgICAgICAgeGF4aXMgID0gbGlzdCh0aXRsZSA9IGxpc3QodGV4dCA9JzxiPlZvcmVzJykpLAogICAgICAgICAgeWF4aXMgPSBsaXN0KHRpdGxlPSdDb3VudCcsIHRleHQ9JzxiPkNvdW50JyksIGJhcm1vZGUgPSAnZ3JvdXAnKQpyX2dyb3VwX2JhcmNoYXJ0CmBgYAojIyMgcGllLWNoYXJ0CmBgYHtyfQpkZl9vcmRlciA9IGRhdGEuZnJhbWUodGFibGUoY2EkcGFya19uYW1lKSkKZGZfb3JkZXIKYGBgCgpgYGB7cn0KcGllX2NoYXJ0ID0gcGxvdF9seSh0eXBlPSdwaWUnLCBsYWJlbHM9ZGZfb3JkZXIkVmFyMSwgdmFsdWVzPWRmX29yZGVyJEZyZXEsIAogICAgICAgICAgICAgICAgICAgIHRleHRpbmZvPSdsYWJlbCtwZXJjZW50JyxpbnNpZGV0ZXh0b3JpZW50YXRpb249J3JhZGlhbCcpICU+JQogICAgICAgICAgICAgICAgICAgIGxheW91dCgKICAgICAgICAgICAgICAgICAgICAgICAgICB0aXRsZT0gbGlzdCh0ZXh0ID0gIjxiPk9yZGVyIERpc3RyaWJ1dGlvbnMiKSwKICAgICAgICAgICAgICAgICAgICAgICAgICBsZWdlbmQgPSBsaXN0KHRpdGxlID0gbGlzdCh0ZXh0PSAnPGI+T3JkZXInKSkpCgpwaWVfY2hhcnQKYGBgCgpgYGB7cn0KaGlzdG9ncmFtX3Bsb3QgPSBwbG90X2x5KGRhdGEgPSBjYSwgeCA9IH4obG9nKHZpc2l0b3JzKSksIG5hbWU9fmNvZGUsdHlwZT0iaGlzdG9ncmFtIikgJT4lCiAgICAgICAgICAgICAgICAgIGxheW91dCgKICAgICAgICAgICAgICAgICAgICAgICAgICB0aXRsZT0gbGlzdCh0ZXh0ID0gIjxiPlRvdGFsIHNsZWVwIHRpbWUgb2YgQW5pbWFscyBiYXNlZCBvbiBWb3JlIiksCiAgICAgICAgICAgICAgICAgICAgICAgICAgbGVnZW5kID0gbGlzdCh0aXRsZSA9IGxpc3QodGV4dD0gJzxiPlZvcmUnKSksCiAgICAgICAgICAgICAgICAgICAgICAgICAgeGF4aXMgID0gbGlzdCh0aXRsZSA9IGxpc3QodGV4dCA9JzxiPlZpc2l0b3JzJykpLAogICAgICAgICAgICAgICAgICAgICAgICAgIHlheGlzICA9IGxpc3QodGl0bGUgPSBsaXN0KHRleHQgPSc8Yj5Db3VudCcpKSkKaGlzdG9ncmFtX3Bsb3QKYGBgCgojIyMgRG9udXQgQ2hhcnQgLyBPcGVuIFBpZS1DaGFydApgYGB7cn0KZGZfdm9yZSA9IGRhdGEuZnJhbWUodGFibGUoY2EkY29kZSkpCmRmX3ZvcmUKYGBgCgpgYGB7cn0KZG9udXRfY2hhcnQgPSBwbG90X2x5KGxhYmVscz1kZl92b3JlJFZhcjEsIHZhbHVlcz1kZl92b3JlJEZyZXEsIAogICAgICAgICAgICAgICAgICAgIHRleHRpbmZvPSdsYWJlbCtwZXJjZW50JykgJT4lCiAgICAgICAgICAgICAgICAgICAgYWRkX3BpZShob2xlID0gMC42KSAlPiUKICAgICAgICAgICAgICAgICAgICBsYXlvdXQoCiAgICAgICAgICAgICAgICAgICAgICAgICAgdGl0bGU9IGxpc3QodGV4dCA9ICI8Yj5PcmRlciBEaXN0cmlidXRpb25zIiksCiAgICAgICAgICAgICAgICAgICAgICAgICAgbGVnZW5kID0gbGlzdCh0aXRsZSA9IGxpc3QodGV4dD0gJzxiPk9yZGVyJykpKQoKZG9udXRfY2hhcnQKYGBgCgpgYGB7cn0KYGBgCgpgYGB7cn0KYGBgCgpgYGB7cn0KYGBgCgpgYGB7cn0KYGBgCgpgYGB7cn0KYGBgCgpgYGB7cn0KYGBgCgpgYGB7cn0KYGBgCgoKYGBge3J9CmBgYAoKCmBgYHtyfQpgYGAK